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ABSTRACT 

We compare the relative merits of AGN selection at X-ray and mid-infrared wavelengths using data from 
moderately deep fields observed by both Chandra and Spitzer. The X-ray-selected AGN sample and associ- 
ated photometric and spectroscopic optical follow-up are drawn from a subset of fields studied as part of the 
Serendipitous Extragalactic X-ray Source Identification (SEXSI) program. Mid-infrared data in these fields are 
derived from targeted and archival Spitzer imaging, and mid-infrared AGN selection is accomplished primarily 
through application of the IRAC color-color AGN 'wedge' selection technique. Nearly all X-ray sources in 
these fields which exhibit clear spectroscopic signatures of AGN activity have mid-infrared colors consistent 
with IRAC AGN selection. These are predominantly the most luminous X-ray sources. X-ray sources that 
lack high-ionization and/or broad lines in their optical spectra are far less likely to be selected as AGN by 
mid-infrared color selection techniques. The fraction of X-ray sources identified as AGN in the mid-infrared 
increases monotonically as the X-ray luminosity increases. Conversely, only 22% of mid-infrared-selected 
AGN are detected at X-ray energies in the moderately deep ((t exp ) sw 100 ks) SEXSI Chandra data. We hy- 
pothesize that IRAC sources with AGN colors that lack X-ray detections are predominantly high-luminosity 
AGN that are obscured and/or lie at high redshift. A stacking analysis of X-ray-undetected sources shows that 
objects in the mid-infrared AGN selection wedge have average X-ray fluxes in the 2-8 keV band three times 
higher than sources that fall outside the wedge. Their X-ray spectra are also harder. The hardness ratio of 
the wedge-selected stack is consistent with moderate intrinsic obscuration, but is not suggestive of a highly 
obscured, Compton-thick source population. It is evident from this comparative study that in order to create 
a complete, unbiased census of supermassive black hole growth and evolution, a combination of sensitive in- 
frared, X-ray and hard X-ray selection is required. We conclude by discussing what samples will be provided 
by upcoming survey missions such as WISE, eROSITA, and NuSTAR. 
Subject headings: galaxies: active - infrared: galaxies - X-rays: galaxies 
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1. INTRODUCTION 

The tight correlation of nuclear black hole mass with the ve- 
locity dispersion and mass of the galactic bulge (e.g., Magor- 
rian et al. 1998; Ferrarese & Merritt 2000; Tremaine et al. 
2002) implies that the growth and evolution of galaxies is 
closely linked to the growth and evolution of the supermas- 
sive black holes which reside in (at least) all massive galaxies 
(e.g., Kauffmann & Haehnelt 2000; Heckman 2008). Indeed, 
recent theoretical work suggests that feedback from active 
galactic nuclei (AGN) plays a dominant role in establishing 
the present-day appearances of galaxies, providing a natural, 
physical explanation for both cosmic downsizing and the pos- 
sibly related bimodality in local galaxy properties (e.g., Scan- 
napieco et al. 2005; Hopkins et al. 2008; Cattaneo et al. 2009). 
However, obtaining an unbiased census of black holes in the 
universe remains challenging, hampering our ability to fully 
probe this connection. 

Most surveys for active galaxies are severely biased towards 
unobscured {type 1) AGN since nuclear emission in such 
sources dominates over host galaxy light at most wavelengths, 

1 Space Radiation Laboratory, Mail Stop 290-17, California Institute of 
Technology, Pasadena, CA 91125 

2 NASA Goddard Space Flight Center, Code 662, Greenbelt, MD 20771 

3 NASA Postdoctoral Program (NPP) Fellow 

4 Columbia University, Department of Astronomy, 550 West 120th Street, 
New York, NY 10027 

5 Steward Observatory, University of Arizona, 933 North Cherry Avenue, 
Tucson, AZ 85721 

6 Jet Propulsion Laboratory, California Institute of Technology, Mail Stop 
169-506, Pasadena, CA 91 109 

7 Quest University Canada, 3200 University Blvd., Squamish, BC 
V8B0N8 CANADA 



making type 1 AGN both more readily identifiable and easier 
to follow up spectroscopically. However, both unified AGN 
models (e.g., Antonucci 1993; Urry & Padovani 1995) and 
the shape of the X-ray background suggest a population of ob- 
scured (type 2) sources which outnumber the type 1 AGN by 
up to a factor of ten (e.g., Comastri et al. 1995; Treister et al. 
2004; Gilli et al. 2007). Determining the ratio of unobscured 
to obscured AGN as a function of luminosity and redshift can 
directly constrain the growth history of supermassive black 
holes, and help to quantify their influence on the evolution of 
their host galaxies (e.g., Ueda et al. 2003; Barger et al. 2005; 
Hasinger et al. 2005; Hopkins et al. 2007). Since different 
search techniques for obscured AGN suffer from different bi- 
ases, the problem is best addressed through multiwavelength 
studies. 

The hard (2-10 keV) X-ray and mid-infrared wavebands 
provide powerful, complementary methods for identifying 
and studying AGN over a wide range of intrinsic obscuration. 
Radiation seen from active galaxies at X-ray energies is pri- 
marily due to direct emission from accretion processes near 
the central supermassive black hole, with higher energy pho- 
tons less susceptible to absorption, while mid-infrared radia- 
tion (rest-frame Ao > 2 pm) is typically dominated by emis- 
sion from dusty obscuring material surrounding the AGN cen- 
tral engine, and is likewise relatively immune to absorption 
(e.g., Rrolik 1999). Hard X-ray and mid-infrared surveys of 
AGN will therefore sample a wide range of intrinsic absorp- 
tion and, in this regard, will be both less biased and more 
complete than surveys in the optical (e.g., Richards et al. 
2006) and soft X-ray (E < 2.4 keV; e.g., Hasinger et al. 1998; 
Schmidt et al. 1998) bands, both of which suffer from strong 
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attenuation by even moderate (A/#> 10 21 cm -2 ) columns of 
dusty, obscuring material. However, AGN samples selected 
in the 2 - 10 keV X-ray and infrared bands will also suffer in- 
completeness. Examples exist of AGN identified from optical 
spectroscopy which remain undetected even in the deepest X- 
ray images yet obtained (e.g., Steidel et al. 2002); similarly, 
there are examples of optically selected AGN not identified in 
infrared surveys (e.g., Stern et al. 2005). Conversely, X-ray 
missions have identified AGN whose optical spectra are de- 
void of AGN signatures (e.g., X-ray bright, optically normal 
galaxies, or "XBONGs"; Comastri et al. 2002). Therefore a 
complete census will necessarily require combining selection 
techniques from several spectral regimes. 

The current generation of X-ray telescopes has dramatically 
advanced our knowledge of the AGN population (see Brandt 
& Hasinger 2005 for a review). The Chandra X-ray Observa- 
tory (Weisskopf et al. 1996) provides a large collecting area, 
a moderate field of view (FOV), and exquisite angular resolu- 
tion (< 1") from 0.5 to 8 keV. These capabilities have allowed 
Chandra extragalactic surveys to efficiently select and opti- 
cally identify large samples of AGN in the 2-8 keV range. 
Observing in this hard X-ray energy band means that even 
sources shrouded by considerable obscuring column densi- 
ties (up to 10 24 cm" 2 ) can be detected at low redshift, and 
sources with even higher column densities can be detected out 
to z > 2. Previous sensitive X-ray telescopes were primarily 
restricted to energies below <~ 2 keV and thus missed many 
of the obscured AGN, a source population that has long been 
theorized to explain the mismatch in spectral shape between 
the 2 - 10 keV X-ray background (T 1.4) and the unob- 
scured active galaxies which dominate the source counts at 
soft X-ray energies (T s=a 1.9; Nandra & Pounds 1994). 

The 2003 launch of the Spitzer Space Telescope (Werner 
et al. 2004) opened a new era in mid-infrared observations, 
providing orders of magnitude improvement in sensitivity in 
the 3.6-160 [im band. The increased sensitivity, combined 
with the large FOV, allows for the first time efficient long- 
wavelength survey capabilities. The primary imaging cam- 
eras on Spitzer are the Infrared Array Camera (IRAC; Fazio 
et al. 2004), providing simultaneous imaging at 3.6, 4.5, 5.8, 
and 8 ^m, and the Multiband Imaging Photometer for Spitzer 
(MIPS; Rieke et al. 2004), providing simultaneous imaging at 
24, 70, and 160 ^m. 

Several methods have been developed to select AGN based 
on their Spitzer colors (Lacy et al. 2004; Stern et al. 2005; 
Alonso-Herrero et al. 2006). These methods exploit the differ- 
ence in the typical spectral energy distribution (SED) of AGN 
compared to 'normal' galaxies. The near-infrared emission 
of both typical galaxies and vigorously star-forming galax- 
ies is primarily produced by a thermal stellar population, re- 
sulting in SEDs peaked near the rest-frame 1.6 fxm "bump." 
This bump, which is caused by the minimum in the opacity 
of the H" ion near 1.6 /im (e.g., John 1988), is a feature of 
almost all stellar populations (e.g., Wright et al. 1994; Simp- 
son & Eisenhardt 1999), whereas AGN-dominated SEDs have 
a non-thermal, roughly power-law shape (for A < 10 fim). 
At longer wavelengths contributions from stellar blackbody 
emission is low, while radiation from the obscuring dust near 
the central engine provides strong, relatively isotropic emis- 
sion which is quite distinct from stellar light (although the 
shape of the AGN SED does somewhat depend on viewing 
angle, e.g., Elitzur 2008). Stern et al. (2005), for example, 
suggested a selection technique exploiting the fact that, for 
AGN, the long-wavelength side of the 1.6 ^m stellar peak 



does not decline; the technique uses an empirically deter- 
mined 'wedge' in IRAC color-color space ([3.6]-[4.5] ver- 
sus [5.8]-[8.0]) that preferentially contains AGN as compared 
to normal galaxies or Galactic stars. Combining a sample 
of 10,000 R < 21.5 spectroscopically identified sources from 
the AGN and Galaxy Evolution Survey (AGES; Kochanek et 
al., in prep.) and mid-infrared observations from the IRAC 
Shallow Survey (Eisenhardt et al. 2004), Stern et al. (2005) 
defined mid-infrared AGN selection criteria which robustly 
identify broad- and narrow-lined AGN, with only 18% sam- 
ple contamination from galaxies (17%) and stars (1%). The 
true sample contamination is likely lower, since many of the 
spectroscopically normal galaxies may harbor active nuclei 
(i.e., are mid-infrared versions of XBONGs). Working from 
the full spectroscopically defined AGES sample, Stern et al. 
(2005) found that the wedge selects 91% of the broad-lined 
AGN, 40% of the narrow-lined AGN, and fewer than 3% of 
the normal galaxies. 

This paper presents an exploration of the relative strengths 
of Chandra and Spitzer as black-hole finders, using a subset of 
six fields from the Serendipitous Extragalactic X-ray Source 
Identification program (SEXSI; Harrison et al. 2003, Eckart 
et al. 2005, and Eckart et al. 2006) for which we have obtained 
mid-infrared coverage with Spitzer. In addition to examining 
the properties of ~ 250 hard X-ray-selected AGN, we extend 
our sensitivity to X-ray emission using a stacking analysis on 
sources with different mid-infrared characteristics. This en- 
ables us to make a comprehensive comparison of the relative 
characteristics of AGN samples selected from their X-ray and 
mid-infrared properties. 

We organize the paper as follows: §2 introduces the SEXSI 
program and presents the complementary mid-infrared ob- 
servations and the resulting source catalog; §3 discusses the 
mid-infrared properties of the 2-10 keV SEXSI sources; 
§4 presents the mean X-ray properties of X-ray-non-detected 
Spitzer sources; §4.3 provides a discussion; and §5 sum- 
marizes our conclusions. Luminosity calculations assume a 
standard cosmology used in our previous work (Eckart et al. 
2006), ft = 0.3, A = 0.7, and H Q = 65 km s" 1 Mpc -1 . 

2. OBSERVATIONS, DATA REDUCTION AND SOURCE CATALOG 

We have expanded the multiwavelength data available for 
six SEXSI fields by obtaining Spitzer imaging observations. 
All six fields have deep Chandra X-ray images, optical imag- 
ing, and extensive, deep optical spectroscopy — all of which 
has been published in Harrison et al. (2003), Eckart et al. 
(2005), and Eckart et al. (2006). In this paper we present 
Spitzer archival and/or targeted observations from 3.6 to 
24 [im from these fields. Below we briefly review the SEXSI 
program, describe the complementary Spitzer sample and data 
reduction, and present the IRAC and MIPS photometry and 
source catalog for all 2-10 keV SEXSI X-ray sources. 

2.1. SEXSI X-Ray Sample and Optical Followup 

The SEXSI survey was designed to obtain a significant sam- 
ple of spectroscopically identified objects in the 2-10 keV 
flux range from a few times 10~ 13 to 10~ 15 erg cm -2 s" 1 , to 
provide a medium-area, medium-depth sample to comple- 
ment both the wide-area, shallow surveys (e.g., HELLAS - La 
Franca et al. 2002; ASCA Large Sky Survey - Akiyama et al. 
2000; ASCA Medium Sensitivity Survey - Akiyama et al. 
2003) and deep, pencil-beam surveys (e.g., CDF-N - Alexan- 
der et al. 2003; Barger et al. 2003; CDF-S - Rosati et al. 
2002; Szokoly et al. 2004). To accomplish this, SEXSI covers 
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more than 2 deg 2 of sky, comprised of twenty-seven archival 
Chandra ACIS images. The resulting source sample consists 
of 1034 hard X-ray-selected sources (Harrison et al. 2003) 
with complete R-bwA imaging (Eckart et al. 2005) and opti- 
cal spectroscopy for 477 (46%) of the sources (Eckart et al. 
2006). The vast majority of these sources are active galaxies 
with L x (2 - 10 keV) > 10 43 erg s _1 . Just over half of these 
AGN have optical spectra that show typical AGN lines which 
we classify as either broad-lined AGN (BLAGN) or narrow- 
lined AGN (NLAGN). Extragalactic sources which exhibit 
optical spectra devoid of the typical high-ionization AGN 
lines are called emission-line galaxies (ELGs) or absorption- 
line galaxies (ALGs); such sources are often referred to as 
XBONGs in the literature. See Eckart et al. (2006) for more 
details on the SEXSI spectroscopic classifications. 

2.2. Spitzer Mid-Infrared Imaging 

We obtained mid-infrared imaging for six SEXSI fields 
through both archival and targeted Spitzer programs (see Ta- 
ble 1). These programs include imaging at 3.6, 4.5, 5.8, and 
8 /im from IRAC and imaging at 24 /im from MIPS. Fields 
were chosen both to optimize use of archival Spitzer data and 
to target SEXSI fields with the most complete spectroscopic 
coverage. Program ID 20808 was designed specifically for 
the SEXSI follow-up, with typical exposure times of 600 s for 
IRAC observations and 1000 s for MIPS observations. Be- 
cause the SEXSI data were drawn from the Chandra public 
archive, the fields are not 'blank' : the target sources in the six 
SEXSI Spitzer fields are all z > 0.18 galaxy clusters, five of 
which have 0.64 < z < 1.3 (see Table 1). 

2.2.1. IRAC Observations and Reductions 

Our IRAC data reduction begins with the basic calibrated 
data (BCD) output from the Spitzer Science Center pipeline 
(ver. 10.0.1). Note that Spitzer BCD images have surface 
brightness units in MJy sr" 1 , where 1 MJy sr" 1 is equal to 
23.5 /iJy arcsec" 2 . Most of the IRAC observations consist of 
six dithered 100 s exposures. For each IRAC channel, the 
BCDs were mosaicked using the MOPEX package (Makovoz 
& Khan 2005) and resampled by a factor of two to give a 
pixel scale of 0.61". The MOPEX outlier (e.g., cosmic ray, 
bad pixel) rejection was optimized for the regions of deepest 
coverage in the centers of the images. Source extraction was 
performed with SExtractor (Bertin & Arnouts 1996) both in 
single-image mode to produce independent source catalogs in 
all four IRAC bands and in dual-image mode using the 3.6 /im 
channel for source detection. All IRAC BCD images have 
astrometry derived from the 2MASS survey and are good to 
better than 0.5". 

In general, we use a 5" diameter aperture for IRAC pho- 
tometry. For the small fraction of the sources (<~ 5%) that are 
spatially resolved and larger than this, we use a 7" diameter 
aperture. We select the apertures to ensure a fixed extrac- 
tion area similar to the isophotal area; this same methodol- 
ogy was adopted by both the Spitzer First Look Survey (FLS; 
Lacy et al. 2005) and the Spitzer High-Redshift Radio Galaxy 
program (SHzRG; Seymour et al. 2007). The aperture cor- 
rections, applied after the conversion from surface brightness 
to flux density, are 1.167, 1.213, 1.237 and 1.466 in the four 
IRAC bands, respectively, for the 5" apertures. The corre- 
sponding corrections for the 7" apertures are smaller (1.100, 
1.130, 1.109 and 1.198, respectively). We derive total un- 
certainties on the resulting IRAC flux densities by adding in 



quadrature the statistical uncertainties (derived from SExtrac- 
tor) to a conservative 10% estimate of the systematic uncer- 
tainties. The latter includes several effects (e.g., Hora et al. 
2008): uncertainties in the IRAC zero point, field-dependent 
color corrections, pixel phase effects (e.g., zero point correc- 
tions that depend on the position of a source relative to the 
nearest pixel center), and variation of the pixel scale across 
the field of view. For nondetections we calculate 5cr upper 
limits. A small number of sources are recovered by SExtrac- 
tor with signal-to-noise ratios below 5 (3 < S/N < 5); exam- 
ination of the images shows that these are valid detections, 
generally near the edge of the mosaic. We flag these sources 
in the catalog and exclude them from the analysis presented 
here. 

2.2.2. MIPS Observations and Reductions 

MIPS observations obtained for this program (e.g., PID 
20808; see Table 1) consist of two cycles of 30 s exposures 
with a 3 x 3 dither pattern at 24 /im. This observing strategy 
provides an average (median) exposure time of 1290 s (877 s) 
across our combined fields of view. To process the MIPS data, 
the BCDs provided by the Spitzer Science Center were first 
median filtered to remove a stripe pattern artifact in the in- 
dividual frames. We then mosaicked the data with MOPEX 
using standard inputs to produce final, stacked images with 
the native platescale of 2.45" pix" 1 . We convert the mosaics 
from surface brightness units to flux density units, and ex- 
tract flux densities with SExtractor using an 1 1 .76" diameter 
extraction aperture with an aperture correction of 1.697. Cal- 
ibration with MIPS is performed with stars using a blackbody 
SED. Fadda et al. (2006) found that the color corrections for 
both starburst galaxies and AGN were small at z < 1 . At z > 1 , 
the color correction for a prototypical luminous AGN (NGC 
1068) remains small. As most SEXSI sources fall within one 
of these categories, we did not apply a color correction to ac- 
count for MIPS sources having non-stellar SEDs. As with the 
IRAC photometry, we derive total uncertainties by adding in 
quadrature the statistical uncertainties (derived from SExtrac- 
tor) to a conservative 10% estimate of the systematic uncer- 
tainties (see, e.g., Engelbracht et al. 2007); for the MIPS data 
the statistical uncertainties are generally the dominant source 
of error. Compared with the IRAC photometry, a larger num- 
ber of MIPS sources are detected with signal-to-noise ratios 
below 5 (3 < S/N < 5). We flag these sources in the cata- 
log and exclude them from our analysis. For the true non- 
detections we calculate 5a upper limits. 

2.3. Source Catalog 

Table 2 presents mid-infrared photometry for 290 hard X- 
ray-selected SEXSI sources. The full machine-readable cat- 
alog is available online. Each of the four IRAC catalogs as 
well as the MIPS catalog were individually matched to the 
SEXSI X-ray source positions 8 using a 2.5" search radius. 
To calculate a false match rate we shifted the X-ray source 
catalog by 1' and matched to the IRAC and MIPS catalogs; 
this entire procedure was repeated six times using different 1' 
shifts. The resulting false match rates are: 10.1% (3.6 /im), 
7.2% (4.5 /im), 3.7% (5.8 /mi), 2.6% (8.0 /im), 1% (24 /zm), 
and < 1% for four-band-detected IRAC sources. For ease of 
reference, columns (4)-(10) of Table 2 present X-ray and op- 
tical source information from previous SEXSI publications, 

8 We use the X-ray positions from Eckart et al. (2005) which have been 
corrected for small Chandra pointing offsets. 
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while the Spitzer mid-infrared data is presented in columns 
(1 1)— (15). Complete X-ray source information is provided in 
Harrison et al. (2003) and optical-followup source catalogs 
are presented in Eckart et al. (2005) and Eckart et al. (2006). 

Column (1) presents the source ID, followed by Aa and 
AS, the difference in right ascension and declination between 
the X-ray position and the Spitzer counterpart. Column (4) 
lists the 2-10 keV flux in units of 10~ 15 erg cm -2 s _1 , while 
column (5) gives the X-ray hardness ratio, HR = (H—S)/(H+ 
S), where H is the hard-band counts (2-10 keV) and S is 
the soft-band counts (0.5-2keV). The optical-counterpart R- 
band magnitude is presented in column (6), followed by the 
optical-spectroscopic classification (7) and the source redshift 
(8). Sources marked 'unid' have no spectroscopic classifica- 
tion or redshift. The observed X-ray luminosity is given in 
column (9) and the N H value determined from the X-ray data 
is found in column (10). When the best-fit Nh measurement 
is zero, we report "<" in the table; see Eckart et al. (2006) for 
details of these calculations. Columns ( 1 1)— ( 14) provide the 
flux density in each of the IRAC bands with associated errors. 
The errors are > 10% of the source flux owing to our conser- 
vative estimate of systematic uncertainties for the Spitzer data. 
The MIPS data are presented in the final column. Entries in 
Columns (1 1)— (15) that are blank indicate that a source was 
not observed in that band. 

3. MID-INFRARED PROPERTIES OF X-RAY SELECTED AGN 

This section explores the mid-infrared properties of the 
sample of 2 - 10 keV SEXSI sources. In particular, we fo- 
cus on whether the X-ray-selected AGN will also be iden- 
tified as AGN via the Stern et al. (2005) IRAC color-color 
AGN selection technique, and how this mid-infrared identi- 
fication depends on various source properties, such as opti- 
cal spectroscopic classification, redshift, X-ray flux and mid- 
infrared flux density (§3.1), X-ray luminosity (§3.2), and R- 
band magnitude (§3.3). Section 3.5 provides a comparison to 
AGN samples presented in the literature: §3.5.1 compares the 
mid-infrared properties of hard X-ray-selected samples, while 
§3.5.2 discusses the identification of high-redshift, obscured 
quasars via 3.6 /rni-24 /im selection criteria. 

The Spitzer data covers 290 hard X-ray SEXSI sources, of 
which 154 have spectroscopic redshifts. However, a signif- 
icant fraction of the X-ray sources do not have coverage in 
all four IRAC bands. 9 For the following analysis we focus 
on the 189 hard X-ray sources (65%) with observations in all 
four IRAC bands. For some analyses, we consider the slightly 
smaller sample of 155 sources (53%) with MIPS 24 ^m ob- 
servations as well. Because of the relative depths of the mid- 
infrared and X-ray data, nearly all of the SEXSI sources have 
infrared counterparts in one or more of the bands. Overall, 
98% of SEXSI sources observed at 3.6 /zm are detected at 
greater than 5 a at that wavelength, and the detection rates for 
the other IRAC bands are all > 90%. Of the four-band IRAC- 
observed SEXSI sources, 57% have optical spectroscopic red- 
shifts, and, of those, 89% are detected in all four IRAC bands. 
The corresponding MIPS detection rate is 54% (68% if < 5a 
detections are included). 

Figure 1 shows an IRAC color-color diagram of all mid- 
infrared sources with > 5a detections in all four IRAC bands 
(small black dots). Larger symbols indicate SEXSI 2 - 10 keV 

9 IRAC has two adjacent fields of view that are imaged in pairs (3.6 and 
5.8 fim; 4.5 and 8.0 fim). The sources in our sample that lack four-band 
coverage typically lie nearer to the edge of the X-ray survey area and were 
imaged by only one IRAC FOV, thus yielding only two-band coverage. 



sources, as detailed in the caption. The dotted lines indicate 
the Stern et al. (2005) AGN selection wedge. We adopt the 
Stern et al. (2005) selection, which has slightly less contami- 
nation than the criteria suggested by Lacy et al. (2004) based 
on [4.5]-[8.0] versus [3.6]-[5.8] color space (e.g., Donley 
et al. 2008). The majority (63%, or 90/142) of the X-ray- 
selected AGN lie in the Stern et al. (2005) selection wedge. 
Considering the spectroscopic sample, the AGN wedge con- 
tains 85% (34/40) of the BLAGN, 64% (9/14) of NLAGN, 
and 34% (13/38) of the ELG, as shown in Figure 1. Since 
there were only two ALGs (and one Galactic star) in the 
SEXSI sample with four-band IRAC observations, these mi- 
nor source populations are not considered further. Includ- 
ing lcr errors on the IRAC flux densities, all but three of 
the BLAGN are consistent with the wedge (one of these is 
the lowest-luminosity BLAGN in the sample, with L2-10 keV = 
4 x 10 42 erg s" 1 ) and all but one of the NLAGN are consistent 
with the wedge selection. In the following sections we explore 
the AGN selection further; in particular, we discuss the dis- 
crepancy between the high fraction of BLAGN and NLAGN 
selected using the mid-infrared techniques as compared to low 
fraction of ELGs. 10 

3.1. Dependence on Observed X-Ray Flux and Mid-Infrared 
Flux Density 

Mid-infrared selection of AGN relies exclusively on ob- 
served mid-infrared colors to distinguish AGN from galax- 
ies. This selection can be compared to X-ray identification 
of AGN for various flux limits in both wavelength regimes. 
Figure 2 explores the wedge membership of our sample for 
flux cuts in the X-ray and mid-infrared. The top panels 
show that X-ray-bright AGN are predominantly found in the 
wedge, while for fainter X-ray sources the wedge member- 
ship fraction decreases. However, the opposite is true when 
considering the mid-infrared data: a significant fraction of 
brighter X-ray-selected AGN are found outside the wedge. 
Roughly half of the SEXSI AGN brighter than / 3 6AJm = 50 fi]y 
([3.6] < 16.86) are outside the wedge, compared to one-third 
for fainter mid-infrared sources. The non- wedge X-ray AGN 
are also generally at low redshift. 

Figure 2 (middle panels) also shows X-ray undetected mid- 
infrared sources. At bright mid-infrared fluxes, most of the 
IRAC sources have colors consistent with low-redshift galax- 
ies (e.g., see galaxy tracks in Figure 1). The X-ray sources 
outside the wedge have similar colors, suggesting that their 
mid-infrared flux is dominated by stellar emission as op- 
posed to AGN activity. Many of these non-wedge-selected 
X-ray sources indeed have optical spectra of ELGs, consis- 
tent with the idea that emission related to star formation out- 
shines the emission from the active nucleus in some parts 
of the spectrum (i.e., in the mid-infrared and optical); how- 
ever, the X-ray luminosities of these sources are too high to 
be attributed wholly to emission related to star formation and 
instead we attribute the X-ray flux to emission from the nu- 
clei. Figure 2 also shows that while a large fraction of the 
brightest mid-infrared IRAC AGN candidates are detected in 
the SEXSI Chandra images, this is not the case at fainter 
mid-infrared magnitudes: there is a large population of mid- 
infrared-selected AGN candidates which lack X-ray detec- 

10 Ten percent of the spectroscopically identified X-ray-selected AGN 
lie at redshifts near that of the target cluster of the Chandra observations. 
Throughout §3 we include these sources in our analysis; however, we note 
that our conclusions remain unchanged if this ten percent of sources is elimi- 
nated from the sample. 
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tions, even in our relatively sensitive ((f eX p) ~ 100 ks) Chan- 
dra data sets. 

Figure 3 shows the mid-infrared and X-ray flux distribu- 
tions of SEXSI sources, and gives an indication of the ef- 
ficiency of mid-infrared AGN selection as a function of the 
imaging depth in both wavelength regimes. The six fields con- 
sidered here cover a relatively small area of the sky, and thus 
there are few bright sources in the sample. However, despite 
these small numbers, our finding that a relatively small pro- 
portion of SEXSI AGN with bright mid-infrared fluxes are se- 
lected by the mid-infrared wedge shows that shallow infrared 
surveys will not be sensitive to the same AGN population as 
X-ray surveys. We consider these results in light of upcoming 
shallow, full-sky survey missions such as eROSITA and the 
Wide-field Infrared Survey Explorer (WISE) in the our con- 
cluding section (§5.2). 

3.2. Dependence on X-ray and Infrared Luminosity 

Figure 4 illustrates the wedge selection for four 2-10 keV 
luminosity ranges, showing that X-ray luminosity strongly af- 
fects the fraction of sources which appear as AGN in the mid- 
infrared. For this analysis, we use the rest-frame 2-10 keV 
X-ray luminosities from Eckart et al. (2006) that were de- 
rived assuming V = 1.5. The entire lowest luminosity sub- 
sample (L x < 10 43 erg s" 1 ) lie either outside of the wedge or 
near the outskirts of the wedge, while the fraction within the 
wedge increases monotonically with X-ray luminosity. At the 
higher X-ray luminosities (L x > 10 43 5 erg s" 1 ), all but two of 
the sources are consistent with the wedge. Of the 36 sources 
with spectroscopic identifications that are not selected by the 
wedge, 25 are classified as ELGs based on their optical spec- 
tra. These sources have a mean redshift of (z) <~ 0.8, average 
infrared luminosity (L4.5) ~ 10 44 erg s" 1 , and average X-ray 
luminosity (L x ) <~ 10 43,2 erg s" 1 . Table 3 tabulates the frac- 
tion of X-ray sources in the wedge for the four plotted X-ray 
luminosity ranges. The X-ray luminosities used for binning 
the sources are the intrinsic, absorption-corrected luminosi- 
ties described in Eckart et al. (2006); these results are not sig- 
nificantly different, however, if observed (e.g., uncorrected) 
luminosities are used instead. 

Overall, it is clear that high X-ray luminosity sources pre- 
dominantly show strong AGN features in their optical spec- 
tra and are wedge-selected, while the low-X-ray-luminosity 
sources generally lack AGN spectroscopic features and are 
not in the wedge. These results show that 2-10 keV X-ray 
surveys are efficient at finding low-luminosity AGN, sources 
missed by mid-infrared selection techniques. This same trend 
is seen when the sample is divided by mid-infrared luminos- 
ity: sources with lower mid-infrared luminosities tend to fall 
out of the wedge. 1 1 This is unsurprising given that, in lumi- 
nous AGN, nuclear emission reprocessed by dust dominates 
the mid-infrared continuum and thus the mid-infrared lumi- 
nosity is strongly tied to the AGN bolometric output. 

3.3. Dependence on R-band Magnitude 

Figure 5 shows this same color-color diagram as in Fig- 
ures 1, 2, and 4 but split at R cut = 21; the top panel shows 

1 1 We do not apply K-corrections when calculating the mid-infrared lumi- 
nosities as they are not well known. We estimate that for quasar-like spec- 
tra the K-correction will not be important, while for galaxy-like spectra it 
will vary a factor of 4-5 between z = 0.5 and z = 1.5. Overall, these K- 
corrections will not change the fact that L x and L/r are correlated, however 
there will be scatter based on the particular SED. 



Spitzer and SEXSI sources with bright optical counterparts 
and the bottom panel shows those with counterparts fainter 
than R cut = 21. Tabulated results comparing the wedge se- 
lection shown in these plots are presented in Table 4. These 
data show a few strong trends. Only 5% ± 1% of optically 
bright four-band IRAC sources lie in the wedge, while 29% 
± 2% of the optically faint four-band IRAC sources do. If 
we consider the subset of these sources that have 2-10 keV 
detections, we see that the 70% ± 16% of the optically bright 
X-ray sources lie in the wedge and 61% ± 7% of the opti- 
cally faint X-ray sources are in the wedge. The X-ray sources 
in general are more likely to fall in the wedge than a typical 
IRAC source, consistent with the notion that the wedge pref- 
erentially identifies AGN (e.g., Gorjian et al. 2008; Ashby 
et al. 2009). An interesting trend evident from this sam- 
ple is the large fraction of R > 21 X-ray undetected wedge 
sources: of the 374 IRAC sources in the wedge with R > 21, 
only 18% ± 2% (68 sources) are X-ray detected. For com- 
parison, the corresponding optically bright (R < 21) sources 
are detected by Chandra nearly two-thirds (63% ± 14%) of 
the time. This trend is consistent with the idea that many of 
the X-ray undetected Spitzer AGN lie at high redshift and/or 
are heavily obscured, resulting in faint optical counterparts 
and X-ray non-detections. The infrared data is deep com- 
pared to the X-ray data, which, in combination with the scat- 
ter in the X-ray to infrared flux ratios, means that sources 
that lie at high redshift may fall below the X-ray detection 
threshold yet have significant detections in the Spitzer data. 
Sources at moderate redshift with unobscured 2-10 keV 
fluxes that would easily be detected in our Chandra data 
will fall below the detection threshold of even the deepest 
Chandra coverage in this study (~ 200 ks, limiting flux of 
«2x 10~ 15 erg cm" 2 s" 1 ) if they are heavily obscured. For 
example, the observed flux of a source at z = 0.4 with an 
unobscured 2-10 keV flux of f x 1 x 10" 14 erg cm" 2 s" 1 
will be reduced to f x = 2x 10" 15 erg cm" 2 s" 1 by an obscuring 
column of 1 x 10 24 cm" 2 , assuming photo-electric absorption 
and T = 1.9. We also note that a significant number of the 
fainter wedge sources may also be normal galaxies (see also, 
e.g., Donley et al. 2008); the typical 600 s IRAC integrations 
considered here reach significantly (« 3x) deeper than the 
shallow data sets from which the original wedge criteria were 
derived. 

The final rows of Table 4 split the population by spec- 
troscopic source type; the columns for Spitzer sources are 
blank here because we do not have spectroscopic followup 
of the non-X-ray-detected Spitzer AGN candidates, though 
follow-up of optically-faint, mid-infrared-selected AGN from 
the IRAC Shallow Survey shows a large fraction of them to 
be type 2 AGN at z > 2 (Stern et al., in prep.). Of the op- 
tically bright ELGs, only one of six sources lies within the 
wedge (and it is at the boundary), while twelve of thirty-two 
R > 21 ELGs lie in the wedge. This trend is different from that 
seen with the BLAGN and NLAGN, where a similar or higher 
fraction of R < 21 sources lie in the wedge. This may indi- 
cate that for optically bright ELGs the starburst component 
of the SED is brighter than the AGN component pushing the 
source away from the AGN-defined wedge, while for some of 
the optically fainter sources the AGN component is strong in 
the mid-infrared leading to wedge selection. This trend is ex- 
plained by higher X-ray luminosity sources tending to fall in 
the wedge: as ELGs are identified at higher redshift they will 
be fainter in R and also have higher average L x due to selection 
effects, and thus preferentially fall in the wedge as compared 
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to low-luminosity, low-redshift sources. Splitting the wedge 
sources by redshift confirms this assertion (Figure 2). Over- 
all, regardless of /^-magnitude, the majority of the ELGs fall 
outside of the wedge. 

3.4. A Brief Summary ofAGN Selection in the X-ray and 
Mid-Infrared 

The results of our comparison of AGN selection techniques 
presented in this section may be summarized as follows: 

• nearly two-thirds (63%) of X-ray-selected AGN lie in 
the mid-infrared selection wedge and nearly all of the 
X-ray-selected AGN that have optical spectra that ex- 
hibit high-ionization emission lines (i.e., BLAGN and 
NLAGN) lie within the mid-infrared selection wedge, 
within photometric uncertainties; 

• on average, the X-ray-selected AGN classified as ELGs 
tend to fall outside of the infrared selection wedge, and 
thus are sources that would be missed by both mid- 
infrared and optical AGN selection techniques; 

• X-ray bright sources tend to fall in the AGN selection 
wedge, while sources with lower X-ray fluxes tend to 
fall outside the wedge; 

• the opposite is true when considering mid-infrared 
fluxes: the X-ray-selected AGN that are faint in the 
mid-infrared tend to fall in the wedge, while those that 
are bright in the mid-infrared tend to be missed by the 
infrared AGN selection technique; 

• the non-wedge X-ray-selected AGN generally lie at low 
redshift; 

• overall, high X-ray luminosity sources fall within the 
mid-infrared selection wedge, while the low-luminosity 
sources fall outside the wedge: the BLAGN and 
NLAGN are predominantly wedge selected, as they 
have higher X-ray luminosities, while the ELGs, which 
have lower X-ray luminosities, tend to fall outside the 
selection wedge; 

• there is a large population of mid-infrared, Spitzer- 
selected AGN that is not detected in the 2 - 10 keV 
X-ray data; 

• when mid-infrared-selected AGN sources are split into 
optically bright (R < 21) and optically faint (R > 21) 
populations, the optically bright sources are detected by 
Chandra the majority of the time, while only w 18% of 
the optically faint sources are identified in the X-ray. 

3.5. Comparison to Previous Work 

3.5.1. Hard X-ray Selection 

We compare the SEXSI Spitzer sample to several other 
Chandra - Spitzer surveys that have a range of depth and 
area. The XMM-Newton Medium Deep Survey (XMDS; 
Pierre et al. 2007) identified 2-10 keV X-ray sources over 
1 deg 2 with mid-infrared coverage from the Spitzer Wide- 
Area Infrared Extragalactic Survey survey (SWIRE; Lons- 
dale et al. 2003). Polletta et al. (2007) and Tajer et al. 
(2007) present a catalog of 122 X-ray-selected AGN with 
fx <; 10" 14 erg cm" 2 s , about a factor five shallower than 
the SEXSI sample considered here. All of these sources have 



at least two-band IRAC detections. Using the Lacy et al. 

(2004) mid-infrared wedge selection technique, they find that 
80-95% of the X-ray sample would be selected on the ba- 
sis of mid-infrared colors. We find that 82% of the XMDS 
AGN are selected by the Stern wedge, much higher than the 
62% of SEXSI sources within the wedge. If a flux cut of 
f x > 10~ 14 erg cm" 2 s" 1 is applied to the SEXSI sources in 
order to match the limit of XMDS, the wedge membership in- 
creases to 74%; when scatter due to photometric errors is con- 
sidered, the wedge membership of the X-ray-brighter SEXSI 
sources is similar to XMDS (see upper panel of Figure 2). 

Barmby et al. (2006) present Spitzer identifications for a 
catalog of ~ 150 X-ray AGN in the Extended Groth Strip. 
The sample includes 138 sources detected in a 200 ks Chan- 
dra exposure. Of these, 67 would be selected on the basis 
of hard-band detection and thus form a comparable sample 
to SEXSI. For the full catalog, Barmby et al. (2006) report 
that only 50% of sources are within the Stern wedge. Con- 
sidering the hard-band selected sample, this percentage in- 
creases to 57%. If a further cut is made to select sources 
with f x >2x 10" 15 erg cm" 2 s" 1 (similar to the typical SEXSI 
limit), the percentage of wedge sources increases to 64%, and 
for a flux cut of /,>lx 10" 14 erg cm" 2 s" 1 (similar to XMDS) 
the percentage reaches 71%. These trends show that various 
surveys agree when similar selection criteria are applied. 

Our results are also consistent with those of Donley et al. 
(2007), who use deeper X-ray data over a smaller area from 
the Ms Chandra Deep Field-North (Alexander et al. 2003) 
and find that a small fraction of low-Z^ sources fit their mid- 
infrared power-law AGN selection criterion, with only one- 
third being within the Stern wedge. The mid-infrared se- 
lection efficiency for their sample increases from ~ 14% at 
io.s-8 keV < 10 42 erg s" 1 to 100% at L .5-8 keV > 10 44 erg s" 1 . 

Combining these results, a consistent picture emerges in 
which hard X-ray surveys uncover a population of low- 
luminosity AGN missed by mid-infrared selection techniques, 
with the efficiency of mid-infrared selection declining in pro- 
portion to the depth of the X-ray data. Optical selection tech- 
niques which rely on the detection of high-ionization narrow 
emission lines in order to identify obscured AGN (e.g., Za- 
kamska et al. 2003) similarly miss this population of low- 
luminosity AGN. 

Lastly, we can compare to the large sample of infrared- 
selected AGN assembled by Hickox et al. (2007) from the 
Bootes field. This sample includes ~ 1500 AGN identified by 
the Stern wedge, with spectroscopic or photometric redshifts 
between 0.7 < z < 3. The shallow X-ray coverage of this field 
results in few X-ray detections for the infrared-selected AGN, 
but X-ray stacking shows that the infrared-selected candidates 
are, on average, consistent with AGN activity. These results 
show that while infrared and X-ray selected samples do not 
identify exactly the same AGN populations, by and large they 
are both effective at constructing samples dominated by AGN. 

We note that broad consistency is found among these vari- 
ous works even though a variety of survey data and selection 
techniques were employed. In particular, the different appli- 
cations of mid-infrared AGN color selection produce similar 
results, albeit with trade-offs in completeness and contami- 
nation. For comparison, we re-analyzed our sample using 
the Lacy et al. (2004) wedge selection technique, and we an- 
alyzed the Tajer et al. (2007) sample using the Stern et al. 

(2005) wedge. We find that the general trends and final con- 
clusions are unchanged, though the numbers and percentages 
do vary. Studies, such as that recently undertaken by Assef 
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et al. (2009), generally find that the Lacy et al. (2004) se- 
lection technique is more complete but is heavily contami- 
nated by low-redshift star-forming galaxies. In comparison, 
the Stern et al. (2005) selection technique is less contami- 
nated by low- and intermediate-redshift galaxies but suffers 
from some incompleteness, notably at z — 4.5 as Ha shifts 
through the IRAC 3.6 fim channel and for sources where the 
AGN emission is dominated by host galaxy stellar emission. 
The power-law AGN selection (Alonso-Herrero et al. 2006; 
Donley et al. 2007) is the least contaminated by non-AGN, 
but only selects a subset of the Stern et al. (2005) and Lacy 
et al. (2004) wedge-selected sources (see, e.g., Donley et al. 
2008). 

3.5.2. High-Redshift, Type-2 Quasar Selection 

Martfnez-Sansigre et al. (2005) use the Spitzer FLS to 
search for type 2 quasars at z > 2 using the following se- 
lection criteria: /24pm > 300 [iJy, f3.(,^ m < 45 ^Jy, and 
350 fiJy < /1.4GHz < 2 mJy. The 24 fim cut is designed 
to define a reliable catalog of luminous quasar candidates, 
while the 3.6 ^m cut eliminates the unobscured quasars and 
the radio cut eliminates starburst galaxies. This effort identi- 
fied 21 high-redshift, type 2 quasar candidates, ten of which 
were spectroscopically confirmed to be narrow-lined AGN at 
1 .4 < z < 4. The other eleven sources exhibit featureless op- 
tical spectra. 

Figure 6 presents our sample of sources detected at both 
24 fim and 3.6 /im, with the upper-left quadrant containing 
the sources that meet the Martfnez-Sansigre et al. (2005) in- 
frared selection criteria. The large symbols indicate SEXSI 
2-10 keV sources, with confirmed z > 2 sources circled in 
black. While our lack of deep radio data precludes applica- 
tion of the radio cut, SEXSI X-ray detections provide another 
effective method to avoid starburst contamination. Of the sub- 
sample of SEXSI sources with MIPS and IRAC coverage, we 
find that all three of the z > 2 confirmed SEXSI NLAGN fall 
in the selection quadrant, while the few confirmed broad-lined 
AGN selected with these criteria all fall along the edges of 
the selection box. The high-redshift, narrow-lined sources ex- 
hibit the high-ionization UV emission lines indicative of an 
active nucleus, and have hard (2-10 keV) X-ray luminosities 
of L x ~ 10 44 5 erg s" 1 . These sources exhibit hard X-ray spec- 
tra, with best-fit Nh estimates of 1 .4-9.4 x 10 23 cm" 2 ; two of 
the three sources do not have significant detections in the 0.5 
- 2 keV band. The sources at z > 2 that exhibit broad AGN 
lines fall outside the quadrant (one sits just at the edge). 

The majority (10 of 17) of the X-ray-selected sources that 
fall in the selection quadrant lack optical spectra. The SEXSI 
optical followup program focused on sources with R < 23 - 
24, and thus many of the fainter 3.6 sources, which tend 
to also be faint in R, lack spectroscopic coverage. However, 
the combined X-ray and optical properties of the unidentified 
sources are consistent with high-luminosity, high-redshift, ob- 
scured AGN. These ten sources are not detected in our /?-band 
imaging or have magnitudes of R > 24 and, if we assume z = 2, 
have high X-ray luminosities (L2-10 keV ^ 10 44 erg s" 1 ). The 
sources are hard in the X-ray: many do not have significant 
detections in the 0.5 - 2 keV band. Using the median hard- 
ness ratio of the sample and assuming z = 2 and an intrinsic 
spectral index of T = 1.9 we find Nh^ 3 x 10 23 cm" 2 . These 
properties are similar to those exhibited by the z > 2 NLAGN 
found by Martfnez-Sansigre et al. (2005), as well as by the 
spectroscopically identified z > 2 NLAGN in our sample. 

We thus confirm the result of Martfnez-Sansigre et al. 



(2005) in a 2 - 10 keV X-ray-selected sample by showing 
that the 3.6 /im-24 /im selection quadrant does indeed select 
high-redshift, highly obscured narrow-lined AGN. Instead of 
using 1 .4 GHz radio selection criterion to eliminate starburst 
contamination we rely on a 2 - 10 keV X-ray detection as an 
effective method to ensure that the sources are active galax- 
ies. All of the spectroscopically confirmed z > 2 narrow-lined 
AGN in the SEXSI sample fall in the selection quadrant, and 
the SEXSI sources in the selection quadrant that lack spec- 
troscopic information have X-ray hardness ratios and /?-band 
characteristics consistent with the notion that they too are 
high-redshift, highly obscured quasars. 

4. X-RAY PROPERTIES OF MID-INFRARED-SELECTED AGN 

A large number of sources in our mid-infrared sample have 
four-band IRAC detections but lack significant Chandra de- 
tections. In this section we present results from an X-ray 
stacking analysis to study the mean X-ray properties of these 
sources, in particular focusing on the subset of these sources 
that fall within the mid-infrared AGN selection wedge. The 
aim is to first understand whether the mean X-ray proper- 
ties of the wedge-selected sources differ from those outside 
the wedge, and whether these differences are consistent with 
the idea that the wedge-selected sources are active galax- 
ies. A second goal is to estimate the intrinsic obscuration of 
the wedge-selected population using the hardness ratio of the 
stacked X-ray signals and assuming various population red- 
shifts. We first describe the stacking methodology, and then 
present the results (§4.2) and a discussion (§4.3). 

4.1. X-ray Stacking Methodology 

The X-ray stacking analysis used tools provided by CIAO 
and was aided greatly by acis_extract 12 (ver. 3.91; 
Broos et al. 2002). We use CIAO (ver. 3.2) and CALDB 
(ver. 3.1) for the analysis presented here. Source lists for X- 
ray spectral extraction were assembled using our SExtractor 
catalogs of four-band IRAC sources (with > 5a detections 
in each band) that fall within the Chandra images, exclud- 
ing sources very near the edge of an ACIS chip. 13 Sources 
that fall within 1 Mpc of the center of each galaxy cluster tar- 
get are eliminated from the stacking to reduce contamination 
from cluster galaxies. 

Spectra were extracted for each source in the resulting cata- 
log (~ 2500 sources). We choose 1.5 keV as the energy at 
which the PSFs are computed and use extraction apertures 
that enclose 80% of the energy (0.8 x PSF). We calculated 
individual auxiliary response files (ARFs) and redistribution 
matrix files (RMFs) for each source, and we extracted a back- 
ground spectrum from a local circular background region that 
includes at least 100 counts, taking care to mask out all de- 
tected X-ray sources (SEXSI sources, including soft-only X- 
ray sources, target point sources, and extended cluster emis- 
sion). We scaled the background spectra based on the ratio 
of total exposure in the source extraction region to that of the 
background region. 

Source counts and scaled background counts were tabulated 
for each source in several energy bands: the standard 0.5-8 
keV (full band), 2 - 8 keV (hard band), and 0.5 - 2 keV (soft 

12 Available at http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html. 

13 These mid-infrared positions were then 'uncorrected': IRAC 3.6 fim 
source positions were shifted to the uncorrected SEXSI X-ray frame by 
adding to each mid-infrared position the average X-ray-to-optical offset cor- 
rection (see §3 and Table 2 of Eckart et al. 2005). 
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band), as well as 2 - 4 keV, 4-6 keV, and 6-8 keV. Esti- 
mates of individual source significance were calculated for 
each band. Following Laird et al. (2005), we define detection 
significance as S/^/B and signal-to-noise ratio as S/VS+B, 
where S and B are the net source counts and scaled back- 
ground counts, respectively. The average exposure per pixel 
(in each energy band) was calculated for each source by aver- 
aging the ARF over the given energy band and then over the 
source extraction cell. In addition to extracting X-ray spectra 
at the position of the IRAC sources, as a control sample we re- 
peated this procedure using a catalog with a 1' offset applied 
to all source positions. 

After calculating individual source statistics we stacked the 
sources. First, any IRAC sources that are within 5" of an 
X-ray source were eliminated from the sample. In addi- 
tion, we eliminated IRAC sources that lie within 5" of 64 
sub-threshold X-ray sources that were initially identified by 
wavelet ect but later eliminated from the main SEXSI cat- 
alogs because they fall below the strict SEXSI P < 10" 6 sig- 
nificance threshold. 14 In field CL0848+44 we include only 
IRAC sources that fall on ACIS chip 3, since the majority of 
the other X-ray sources are at very large off-axis angles. 

The total signal in each band was calculated by summing 
the net background-subtracted counts from each source in the 
sample. To calculate an average flux we divided the summed 
signal by the sum of the average exposure per pixel in each ex- 
traction cell to convert from counts to ph cm" 2 s" 1 . To convert 
from ph cm" 2 s" 1 to erg cm" 2 s" 1 , we assumed a power-law 
spectrum with photon index T = 1 .5. 

4.2. Stacking Results 

We stacked three basic samples: all X-ray undetected 
sources (1298 sources), the wedge-selected sources (226 
sources), and those that fall outside of the wedge (1072 
sources). Table 5 presents a summary of the stacking re- 
sults, where N is the number of sources in the particular sam- 
ple. The "Cts" column gives the total number of background- 
subtracted net counts in the sample. The "Sig" column 
gives an estimate of the significance of the stacked counts 
(Stot/V^tot, see §4-1) and f x is the average X-ray flux for that 
band in units of 10~ 17 erg cm" 2 s" 1 . The full sample of all 
X-ray undetected, four-band IRAC sources showed a signifi- 
cant detection in the stacked Chandra images. The stacks of 
the sources created by the shifted source catalog showed no 
significant detection in any band. 

Stacking only the wedge-selected sources, we find signifi- 
cant detections in all three X-ray bands. The wedge-selected 
sample is, on average, brighter in the X-ray: the average 
0.5 - 8 keV X-ray flux from the wedge-selected IRAC sam- 
ple is ~ 2 times larger than that of the full sample of four- 
band IRAC sources. This trend is most significant in the hard 
band, where the average flux of the wedge-selected sample is 
over three times that of the sample that lies outside the wedge. 
These findings are consistent with the idea that AGN prefer- 
entially lie inside the wedge. 

To allow for errors in the IRAC photometry (e.g., see Fig- 
ure 1), we also stacked based on an extended wedge defined 

14 In Harrison et al. (2003) we used wavdetect for initial source iden- 
tification with the sigthresh parameter set to 10~ 5 to produce a complete 
catalog. In a subsequent step we tested the significance of each source in- 
dividually and eliminated sources from the source catalog with a nominal 
chance occurrence probability greater than l(r 6 . Sources that fall between 
the two thresholds were not included in the SEXSI X-ray sample. 



by: 

([5.8]- [8.0]) > 0.5 
A ([3.6]-[4.5])>0.2-([5.8]-[8.0]) + 0.078 (1) 
A ([3.6]-[4.5])>2.5-([5.8]-[8.0])-3.77, 

where A is the logical AND operator. This extended wedge 
is constructed by taking the original Stern et al. (2005) wedge 
and adding a swath of width 0. 1 mag at each edge. The stack- 
ing results for the samples inside and outside the extended 
wedge are also presented in Table 5; the extended wedge 
swath includes an additional 150 sources. Removing these 
150 sources from the 1072 sources not selected by the origi- 
nal wedge removes a third of the stacked counts in the hard- 
band, leaving a source significance of only 3.9 in the hard- 
band sample of sources outside of the extended wedge. This 
exercise demonstrates that the hard X-ray emission, which is 
more easily produced by AGN than by star formation, comes 
primarily from sources in or just outside the wedge, at least at 
our detection sensitivity. 

Comparing the stacked fluxes of wedge-selected sources to 
those outside the wedge we find that the mid-infrared-selected 
AGN are about three times brighter in the hard band than 
those outside the wedge, which is the same ratio we find if 
we calculate the average flux for SEXSI 2-10 keV-detected 
sources inside and outside the wedge. Although the wedge- 
selected sources are harder than those outside the wedge, the 
column density implied by the average hardness ratio is not 
extreme. Table 6 provides estimated Nh values based on the 
hardness ratio of the stacked spectrum of the wedge-selected 
sample for assumed population redshifts ranging from 0.0 to 
3.0. The Nh values are calculated assuming an underlying 
power-law spectrum with T = 1 .9 and a Galactic column den- 
sity of 2 x 10 20 cm" 2 ; the associated uncertainties are esti- 
mated by propagating the Poisson sampling error in the hard- 
ness ratio. These Nh estimates range from <~ 1.3 x 10 22 to 
4.8 x 10 23 cm" 2 , suggesting that the sources are intrinsically 
obscured, but are not, on average, Compton thick. 

4.3. Discussion of Stacking Results 

The hardness ratio of the wedge-selected stacked spectrum 
is consistent with moderate intrinsic obscuration, but is not 
suggestive of a highly obscured, Compton-thick source pop- 
ulation. It is possible that the stacking procedure is miss- 
ing a population of sources, such that the stacked counts are 
dominated by a subset of the less obscured wedge-selected 
sources. Another possibility is that the assumption of an 
intrinsic T = 1 .9 power-law component plus photoelectric ab- 
sorption is too simple, so that the estimated Nh values are not 
representative of the sample. For example, it is possible that 
the sources within the sample have a range of power-law in- 
dices and/or absorbing columns. In addition, the soft X-ray 
emission may originate in a different location than the hard, 
power-law component, skewing the Nh estimates to lower val- 
ues. In reality, it is likely that some combination of these and 
other effects is operating. 

Hickox et al. (2007) use an IRAC color-selection technique 
to identify a large sample of AGN from the IRAC Shallow 

15 We estimate that a lower limit of 30% of the sources must be con- 
tributing to the stacked signal: if it were fewer than 30% then individual 
source fluxes would fall above the Chandra detection threshold in order to 
reproduce the average stacked flux in the 2 — 8 keV band. This provides an 
upper limit of 70% of the sources that are not contributing significantly to the 
stacked signal. 
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Survey, and an infrared/optical selection technique to sep- 
arate the sources into unobscured (IRAGN1) and obscured 
(IRAGN2) sub-samples. They perform a similar stacking 
analysis of complementary X-ray data from the shallow (5 ks) 
Chandra XBootes survey and find that the IRAGN1 are con- 
sistent with no obscuration and the IRAGN2 are consistent 
with Nh~ 3 x 10 22 cm" 2 , independent of photometric red- 
shift from Zphot ~ 0.7-2.5. The implied obscuring column 
density found using this shallow, wide-area data set is con- 
sistent with that of the SEXSI stacking analysis assuming 
(zstack) ^1- If the majority of the stacked emission in the 
SEXSI wedge-selected sample comes from a higher redshift 
population, then the implied Nh is higher than that found 
in the XBootes analysis. The average X-ray fluxes of X- 
ray undetected infrared AGN from Hickox et al. (2007) are 
fi-i kev ~ 2 x 10~ 15 erg cm" 2 s" 1 , which is just at the hard- 
band SEXSI X-ray detection threshold, and more than an or- 
der of magnitude larger than the average stacked hard-band 
fluxes from the SEXSI analysis. 

Several groups have explored the X-ray properties of small 
samples of optically identified type 2 quasars from the Sloan 
Digital Sky Survey using Chandra and XMM-Newton (e.g., 
Vignali et al. 2006; Ptak et al. 2006). These studies have 
found that approximately half of the type 2 quasars are de- 
tected in the w 5 - 40 ks X-ray images, and those that are 
detected exhibit moderate obscuring column densities (~ 
10 22 - 10 23 cm ). The sources that are undetected in the 
X-ray are likely highly obscured such that little X-ray radi- 
ation below 10 keV escapes. The high bolometric luminosi- 
ties of these sources, as evidenced by both their mid-infrared 
properties and their optical spectra, imply highly obscured, 
high-luminosity AGN that will require a sensitive high-energy 
(> 10 keV) X-ray mission before they are detected in the X- 
rays. 

5. SUMMARY AND CONCLUSIONS 

We have presented an exploration of the relative efficiency 
of selecting AGN in the X-ray with Chandra and in the 
mid-infrared with Spitzer. In this section we summarize our 
comparison of the selection techniques and make predictions 
about the AGN samples that will be selected with future X-ray 
and mid-infrared surveys. 

5.1. Summary 

Two-thirds of the X-ray-selected AGN are also selected by 
the Stern et al. (2005) Spitzer AGN selection criteria. Nearly 
all of the SEXSI X-ray sources that are classified as either 
broad-line or narrow-line AGN (i.e., those that exhibit high- 
ionization optical/UV emission lines in their optical spectra) 
are consistent with wedge selection within photometric uncer- 
tainties. 

A large fraction of the low luminosity (L2-10 kev < 
lO^ergs" 1 ), hard X-ray (2 - 10 keV) AGN identified by 
Chandra will be missed by Spitzer selection techniques. 
These low X-ray luminosity sources tend to lack the high- 
ionization emission lines that allow AGN identification via 
optical spectroscopy; thus, the sources that are missed by the 
Spitzer color-color AGN selection will also remain unidenti- 
fied in optical spectroscopic surveys. These sources tend to lie 
at lower redshift and have higher fluxes in the infrared, on av- 
erage. This suggests that many of these non-wedge-selected 
X-ray sources have mid-infrared fluxes that are dominated by 
stellar (starburst) emission as opposed to AGN activity; in- 
deed, most have optical spectra of ELGs, consistent with the 



conclusion that emission related to star formation outshines 
the emission from the active nucleus in some parts of the spec- 
trum (i.e., in the mid-infrared and optical), while the X-ray 
luminosities of these sources are dominated by X-ray emis- 
sion from the active nuclei. In addition, we find that while a 
fair fraction (<~ 63%) of the optically bright (R < 21) wedge- 
selected Spitzer AGN are also identified by Chandra, over 80 
percent of the optically faint wedge sources are not detected 
in the X-ray. It is likely that these fainter sources are AGN 
that are heavily obscured and/or lie at high redshift. We also 
find that X-ray-selected high-redshift, type-2 quasars are also 
selected via the 3.6 /xm-24 selection criteria proposed by 
Martmez-Sansigre et al. (2005); the radio selection criteria 
used by those authors to eliminate starburst contaminants is 
not necessary when the sources have X-ray properties indica- 
tive of AGN activity. 

To explore the mean X-ray properties of the X-ray unde- 
tected, IRAC-selected AGN candidates, we have stacked the 
Chandra data extracted from the positions of IRAC sources. 
The stacked wedge sources show significant X-ray signals in 
the full, soft, and hard X-ray bands. The average hardness 
ratio of the stacked spectrum of sources inside the wedge is 
higher than that of sources outside of the wedge, and the aver- 
age 2-8 keV flux is three times larger. The hardness ratio of 
the wedge-selected stacked spectrum is consistent with mod- 
erate intrinsic obscuration 10 22 - 10 23,7 cm" 2 ), but is not 
suggestive of a highly obscured, Compton-thick source popu- 
lation. 

5.2. Conclusions 

The results of our study illustrate the challenges of identi- 
fying obscured and low-luminosity AGN using any one tech- 
nique. A complete understanding of AGN activity and black 
hole growth will require multiwavelength data sets: low- 
luminosity AGN in galaxies with active star formation cannot 
be selected using mid-infrared color techniques, but are ef- 
fectively found at low redshift with the current-generation of 
X-ray telescopes; obscured high-redshift objects are beyond 
the sensitivity threshold of Chandra and XMM, but they are 
effectively identified in the mid-infrared if the AGN luminos- 
ity is at the high end of the luminosity distribution. 

The results from current X-ray and mid-infrared satellites 
enable us to estimate the progress possible with upcoming 
space missions that plan to conduct all-sky or pencil beam 
surveys. The Wide-field Infrared Survey Explorer (WISE; 
Mainzer et al. 2005), scheduled for launch at the end of 2009, 
will provide full-sky mid-infrared imaging reaching a depth of 
120 /xJy at 3.3 /xm (similar to IRAC channel 1) and 160 [iJy at 
4.7 /im (similar to IRAC channel 2). At these shallow depths, 
almost every source redder than [3.3]-[4.7] « 0.5 will be an 
AGN; high-redshift galaxies are too faint to make this flux cut, 
and the coolest brown dwarfs are too rare. The primary con- 
taminant will be actively star-forming galaxies at redshifts of 
a few tenths, which should be readily identifiable from their 
optical morphologies and colors in relatively shallow imag- 
ing such as is available from the Sloan Digital Sky Survey. 
By applying a flux cut at the WISE depths to the Bootes 
IRAC survey (Eisenhardt et al. 2004; Ashby et al. 2009), 
we predict that WISE should detect approximately 85 AGN 
candidates per square degree, of which m 15% will be low- 
redshift star-forming galaxy contaminants. WISE will there- 
fore provide an unprecedented sample of high-luminosity ob- 
scured AGN across a range of redshifts up to z ~ 3 identified 
based on mid-infrared color selection alone. WISE will also 
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detect numerous other low-redshift (z < 1), low-luminosity 
(£2-10 keV ^ 10 43 ergs _1 ) active galaxies, but the significant 
number of open symbols above the WISE threshold shown in 
Figure 3 shows that mid-infrared colors alone are insufficient 
to identify them as AGN. 

The extended Rontgen Survey with an Imaging Telescope 
Array (eROSITA; Predehl et al. 2006) X-ray mission, planned 
for launch in 2011, will provide full-sky hard X-ray images 
to a depth of /2-10 kev = 1-5 x 10" 13 erg cm" 2 s" 1 . Using the 
sample of XBootes sources (Kenter et al. 2005) with fluxes 
above the eROSITA hard-band sensitivity limit, we predict that 
eROSITA should detect w 4 hard X-ray sources per square de- 
gree. The mission will detect many more AGN in the soft 
X-ray band because the effective area is much larger at en- 
ergies below 2 keV, reaching soft-band depths of /0.5-2 kev = 
9 x 10" 15 erg cm" 2 s" 1 . Scaling from the XBootes 0.5 - 2 keV 
sample we predict that eROSITA will detect « 120 sources per 
square degree in the soft band. Figure 3 shows that eROSITA 
will identify a large population of AGN below the WISE 
threshold, most of which we project will be at z > 1. The 
majority of these sources will be unobscured BLAGN, but 
a significant population will be narrow-line, obscured AGN 
(NLAGN and ELG). However, eROSITA does not reach suffi- 
cient depth to detect many of the low-luminosity, low-redshift 
AGN detected, but not identified as such, by WISE. 

Clearly WISE and eROSITA combined will miss (or mis- 
identify) low-luminosity obscured AGN at all redshifts, most 
of which will also not be found in optical and soft X-ray 
surveys. At low redshift only a sensitive X-ray mission ex- 
tending to E x > 10 keV will effectively uncover this popu- 
lation, which will be visible only in reprocessed light, if at 
all, at E x < 10 keV. The Nuclear Spectroscopic Telescope Ar- 
ray (NuSTAR; Harrison et al. 2005) 16 , scheduled for launch 
in 2011, will make major advances in studying faint extra- 
galactic sources at energies above 10 keV. NuSTAR will be 
the first focusing, high-energy X-ray satellite and is designed 
to make targeted observations of the hard X-ray sky from 
6-79 keV. NuSTAR, which has a field-of-view comparable 
to Chandra, plans to cover the Bootes and Great Observato- 
ries Origins Deep Survey (GOODS) fields, and should resolve 
> 50% of the 30 keV X-ray background, reaching 10-40 keV 
flux limits of < 10" 14 erg cm" 2 s" 1 . Many of the resolved 
sources are expected to be low-luminosity, nearby AGN not 
identified as such in current mid-infrared or 2 - 10 keV X- 
ray surveys, and some NuSTAR sources may be AGN with 
hard spectra owing to low-efficiency accretion. For study- 
ing the high-redshift, low-luminosity population, the Interna- 
tional X-ray Observatory (IXO) will reach 2-10 keV flux 
limits of <~ 5 x 10" 17 erg cm" 2 s" 1 in a megasecond exposure 
with its wide-field imager. It could survey up to 10 deg 2 to 
a depth of 10" 16 erg cm" 2 s" 1 and 0.25 deg to a depth of 
10" 17 erg cm" 2 s" 1 , uncovering the largest population of high- 
redshift AGN yet detected. 

Predictions for the AGN populations that will be uncovered 
by these future missions must rely on careful studies of the 
AGN identified by existing X-ray and infrared surveys. As il- 
lustrated in this paper, an understanding of the AGN selection 
techniques of the various missions, and their relative strengths 
and weaknesses, is crucial as datasets are synthesized to ob- 
tain an understanding of black hole growth and evolution in 
the universe. 

Finally, we attempt to put the results of our study in the 

16 http : / /www .nustar.caltech.edu. 



context of understanding the cosmic history of black hole 
growth and AGN accretion modes. The primary conclusions 
of this paper are that heavily obscured luminous AGN are of- 
ten missed by X-ray selection, while low-luminosity AGN 
are often missed by mid-infrared selection. For the high- 
luminosity sources, various work has shown that, at least 
for the optically-selected unobscured sources, most are emit- 
ting at close to their Eddington limit (e.g., Kollmeier et al. 
2006). Multiwavelength observations are consistent with the 
heavily-obscured high-luminosity sources having similar in- 
trinsic SEDs, simply with their UV, optical, and X-ray emis- 
sion suffering from absorption. 

As for the low-luminosity AGN, there are two possibili- 
ties: they could be lower mass black holes again emitting 
at close to their Eddington limits, or they could be average- 
sized black holes, with masses similar to the high-luminosity 
sources, but simply emitting at lower Eddington ratios. Us- 
ing a new set of empirical, low-resolution SED templates for 
AGN and galaxies, Assef et al. (2009) shows that the likeli- 
hood of a source to be selected as an AGN using the Stern 
et al. (2005) mid-infrared color criteria is primarily a func- 
tion of how strong the AGN is relative to the host galaxy — 
e.g., the Eddington ratio assuming that, to first order, the host 
galaxy mid-infrared emission scales with the galaxy total stel- 
lar mass which scales with the nuclear supermassive black 
hole mass. Thus, the fact that we find few low-luminosity 
sources using the mid-infrared criteria suggests that the low- 
luminosity sources are not emitting at close to their Edding- 
ton ratio. This is consistent with the results of Babic et al. 
(2007) who find that low-luminosity X-ray sources in the 
Chandra Deep Field-South have a wide range of Eddington 
ratios, 10" 5 < log(Lboi/^Edd) ^ L One caveat on this interpre- 
tation is that it assumes that the AGN SED is independent of 
both Eddington ratio and luminosity — the former of which 
has recently been questioned by Vasudevan & Fabian (2007). 
Future work, in particular, higher energy observations with 
NuSTAR, will quantify the extent to which the intrinsic AGN 
SEDs depend on both luminosity and Eddington ratio. 

This work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Lab- 
oratory, California Institute of the Technology under con- 
tract with the National Aeronautics and Space Administra- 
tion (NASA). Support for this work was provided by NASA 
through award number 1314516 issued by JPL/Caltech. The 
authors thank the anonymous referee for his/her careful read 
and unusually diligent comments as well as Lewis Kotredes 
for assistance with Spitzer data reduction. MEE acknowl- 
edges support from the NASA Postdoctoral Program. 

Facilities: Spitzer, CXO, Keck:I, and Keck:II. 
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TABLE 1 

Spitzer OBSERVATIONS OF SEXSI FIELDS. 



Cluster Field 


z 


Chandra 
Exp. Time 3 


Spitzer 
Instrument 


PID b 


notes 


f -• \ f\rt A f> . 4 A 

CL0848+44 


1.27 


186 ks 


IRAC 


00017 










IRAC 


00064 










IRAC 


20694 


wider FOV 








MIPS 


00083 










MIPS 


20694 


wider FOV 


RXI 0910+54 


1.11 


171 ks 


IRAC 


00017 










IRAC 


20808 


wider FOV 








MIPS 


00083 










MIPS 


20808 


wider FOV 


RXJ1317+29 


0.81 


111 ks 


IRAC 


00017 




BD1338+29 


0.64 


38 ks 


IRAC 


20808 










MIPS 


20808 




RXJ1716+67 


0.81 


52 ks 


IRAC 


00017 










IRAC 


20808 


wider FOV 








MIPS 


20808 




RXJ2247+03 


0.18 


49 ks 


IRAC 


20808 










MIPS 


20808 





"This column provides the Chandra exposure times after removal of back- 
ground flares. 

'"Spitzer program identifier (PID) 00017 — distant X-ray galaxy clusters 
(P.I. Fazio; GTO); PID 00064 — combined program (P.I. Fazio; GTO); 
PID 00083 — use of massive clusters as cosmological lenses (P.I. Rieke; 
GTO); PID 20694 — IRAC and MIPS mapping of galaxy populations in 
a supercluster at z = 1.27 (P.I. Stanford; G02); PID 20808 — SEXSI (P.I. 
Stern; G02). 



TABLE 2 
SEXSI-Spitzer Catalog. 





A™ 




■J2-10 


HI? 
11 1\ 


A 


class 


z 


\r,n(J \ 
YOg(L x ) 




J 3.6 


J 4.5 


J5.8 


78.0 






I J 


L J 


r i n— 1 5 i 
[10 J 




[mag] 












[fjjy] 


[^Jy] 


[^Jy] 


[/Jy] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(') 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


J084809. 8+444901 


-0.7 


0.1 


35.50 


-0.26 


20.73 


unid 








82. 5± 8.3 




133.9± 13.6 






J08481 1.7+445302 


-1.7 


0.1 


8.01 


0.44 


23.46 


unid 








49. 8± 5.0 


48. 6± 4.9 


48.0± 5.3 


55. 7± 6.1 




J084812.4+445657 


0.3 


0.5 


15.40 


-0.32 


19.45 


unid 










210.2±21.1 




584.4± 58.6 




J0848 18.4+444844 


-0.9 


-0.0 


9.18 


-0.40 


20.66 


ELG 


0.405 


42.7 


< 


65. 7± 6.6 




47.4± 4.9 






J084820. 8+445648 


1.8 


-1.5 


4.55 


-0.15 


>24.40 


unid 










<8.5 




<18.0 




J084822.2+445223 


-2.0 


-0.7 


3.54 


-0.04 


23.32 


BLAGN 


2.187 


43.9 


20.8 










< 76.0 


J084822.2+445627 


-1.5 


-1.2 


10.90 


-0.26 




unid 










346.1 ±34.7 




173.1± 17.9 




J084824.6+445355 


-0.8 


0.3 


2.61 


0.03 


22.21 


ELG 


0.747 


42.7 


22.4 










711. 8± 101.1 


J084825. 1+444808 


-1.7 


-0.2 


7.61 


-0.55 


20.73 


BLAGN 


1.320 


43.7 


< 


84.5± 8.5 




104.5± 10.6 






J084827.2+445433 


-0.6 


-1.0 


9.98 


-0.40 


20.33 


BLAGN 


0.899 


43.5 


< 










1228.5± 155.5 


J084827.4+445604 


2.1 


-1.6 


4.56 


0.31 


23.67 


ELG 


1.528 


43.7 


22.7 










564.4± 86.1 


J084830.2+445604 


0.0 


0.0 


2.55 


-0.18 


>24.40 


cont 
















<65.9 


J08483 1.6+445442 


-0.3 


0.7 


2.44 


-0.15 


25.42 


ELG 


1.267 


43.2 


22.0 


11.8± 1.2 


11.5± 1.2 


10.6± 1.6 


8.5± 1.3 


< 90.6 


J084832.7+445711 


0.0 


0.0 


2.62 


0.57 


22.01 


ELG 


0.749 


42.7 


23.0 










< 87.1 


J084834.4+444937 


1.2 


0.0 


3.02 


0.15 


>24.40 


unid 








9.5± 1.0 




12.0± 1.7 




< 104.2 



Note. — 



Columns (4)-(10) present X-ray and optical source information from previous SEXSI publications ( Harrison el al. 2003 : Eekarl el al. 2005: Eckarl el al. 2006). 
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TABLE 3 

Wedge selection as a function of X-ray luminosity. 



Sample 


log(imin) log(Lmax) 


Wedge 


Total 


% 


All classes 


43.0 


5 


14 


36% ± 16% 


BLAGN 




1 


2 


50% 


NLAGN 







1 


0% 


ELG 




4 


11 


36% 


— ATI — 1 

All classes 


XTTt TT~c 

43.0 43.5 


13 


26 


50% ± 14% 


BLAGN 




4 


5 


80% 


NLAGN 




2 


2 


100% 


ELG 




7 


19 


37% 


All classes 


43.5 44.0 


14 


23 


61% ± 16% 


BLAGN 




11 


14 


79% 


NLAGN 




1 


2 


50% 


ELG 




2 


7 


29% 


All classes 


44.0 


23 


28 


82% ± 17% 


BLAGN 




17 


18 


94% 


NLAGN 




6 


9 


67% 


ELG 







1 


0% 


Note. — 


The percentage errors represent the Poisson sampling un- 



certainties. 



TABLE 4 

Wedge selection as a function of optical magnitude. 





Bright: R < 21 




Faint: R > 21 








All 4-band IRAC 


SEXSI X-ray 


All 4-band IRAC 


SEXSI X-ray 




Sample 


Wedge Total % Wedge 


Total % 


Wedge Total % Wedge 


Total 


% 


All 


30 638 5% ± 1% 19 


27 70% ± 16% 


374 1298 29% ±2% 68 


111 61 


% ± 7% 


BLAGN 


13 


13 100% 


20 


26 


77% 


NLAGN 


2 


3 67% 


7 


11 


64% 


ELG 


1 


6 17% 


12 


32 


38% 


No ID 


3 


5 60% 


29 


42 


69% 



NOTE. — The percentage errors represent the Poisson sampling uncertainties. 



TABLE 5 

Average properties of X-ray undetected IRAC sources. 



Sample 


N 


0.5 


-8keV 


0.5 - 2 keV 


2 


-8keV 




Cts 


Sig. 


fx 


Cts 


Sig. 


fx 


Cts 


Sig. 


fx 










[io- 17 ] 






[io- 17 ] 






[io- 17 ] 


All 


1298 


699.2 


18.7 


10.4 


427.3 


19.8 


2.1 


271.9 


8.9 


9.2 


Inside Wedge 


226 


243.5 


14.5 


19.1 


123.8 


12.8 


3.3 


119.6 


8.7 


21.4 


Outside Wedge 


1072 


455.7 


13.6 


8.4 


303.4 


15.7 


1.9 


152.3 


5.6 


6.4 


Inside Extended Wedge 


376 


367.2 


16.9 


17.1 


192.6 


15.5 


3.0 


174.6 


9.8 


18.5 


Outside Extended Wedge 


922 


332.0 


10.9 


7.2 


234.7 


13.3 


1.7 


97.3 


3.9 


4.8 



NOTE. — N indicates the number of sources in the particular sample, Cts provides the number of net background- 
subtracted counts in each stacked sample, Sig. provides an estimate of the source significance (S/VB), and f x is the 
average X-ray flux for the given band in units of 10~ 17 erg cirr 2 s -1 . See §4.1. 
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TABLE 6 
Wedge- 
selected 
stacked 
spectrum: 

IMPLIED Nh 



z N H 

[10 22 cm" 2 ] 

0.0 1.3 ±0.2 

0.5 3.5 ±0.6 

1.0 7.4 ±1.4 

1.5 14±3 

2.0 24 ±5 

3.0 48 ±9 
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FIG. 1 . — Color-color diagram of sources detected at > 5c in all four IRAC bands (dots). Sources with a 2 - 10 keV X-ray detection have a larger symbol 
overlaid, coded based on the optical spectroscopic classification (see inset). The dotted lines demarcate the AGN selection wedge introduced by Stern et al. 
(2005). The < z < 3 color tracks for two non-evolving galaxy templates from Devriendt et al. (1999) are illustrated; the large filled circles indicate z = 0. A 
starburst galaxy is illustrated with the track of M82 (dashed line) and NGC 4429, an SO/Sa galaxy with a star-formation rate ~ 4000 X lower, is indicated with 
the solid line. Median errors are illustrated at the lower left. 
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FIG. 2. — Mid-infrared wedge selection for samples divided by hard X-ray flux (upper panels; units of erg cm -2 s _1 ), 3.6 /im flux density (middle panels; 
Vega magnitudes) and redshift (lower panels). The upper left panel shows a median error bar for each color. Wedge selection is much more efficient at higher 
X-ray fluxes, as the brighter X-ray sources are predominantly BLAGN. However, wedge selection is not more efficient at higher mid-infrared flux densities; a 
significant number of SEXSI ELG AGN outside of the wedge are brighter than [3.6] = 17 (44 fi]y). The background dots in the middle panel show all four-band 
IRAC sources in the SEXSVSpitzer data; comparison to Figure 1 shows that the brighter sources have IRAC colors similar to starburst galaxies. The lower panel 
shows that the sources missed by the wedge are mainly at low redshift. Taken together, these results imply that shallow mid-infrared surveys would detect a large 
number of AGN that have galaxy colors in the infrared, do not show strong AGN lines in their optical spectra, and are X-ray faint. These sources can only be 
readily identified with deep X-ray surveys. 
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FIG. 3 . — Mid-infrared flux density vs. hard X-ray flux for SEXSI sources. Filled symbols indicate sources within the Stern wedge. The histogram on the right 
shows the 3.6 fim flux distribution for SEXSI AGN (light gray) and all IRAC wedge sources from the SEXSI Spitzer fields (black). The dotted lines illustrate 
the depths that will be reached by two upcoming full-sky survey satellites: the Wide-field Infrared Survey Explorer (WISE; horizontal line) and eROSITA (vertical 
lines). The line labeled "eROSITA (hard)" provides the expected 2-10 keV sensitivity while the "eROSITA (soft)" line gives an indication of the 0.5 - 2 keV 
sensitivity, scaled to the 2-10 keV-band using T = 1.9. The eROSITA effective area is much larger below 2 keV, thus the mission will detect many more AGN in 
the soft band. Some SEXSI AGN brighter than the "eROSITA (soft)" limit but with spectra harder than T = 1.9 will in fact not be detected by eROSITA. 
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FIG. 4. — IRAC color-color diagram of SEXSI 2-10 keV sources, split by absorption-corrected intrinsic X-ray luminosity. As the X-ray luminosity grows, a 
higher fraction of sources fall solidly into the Stern et al. (2005) mid-infrared AGN selection wedge. 
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FIG. 5. — IRAC color-color diagram of sources with optical counterparts brighter (top panel) and fainter (bottom panel) than R = 21. The galaxy tracks are the 
same as in Figure 1 . 



20 



3.5 - 



=1 



V) 

c 

Q 

X 
3 



E 

CM 



3.0 - 



2.5 - 



? 2.0 



1.5 



□ BLAGN 
ONLAGN 
AELG 
no ID 



O z > 2 



A- 

















■ '. ' o ' 






■ °A, 






6 




<g 








A 








o 



I a P. 

: A 

<» S 

Q □ A 



■ 

□ 



A 

o 



■ n _ z 

0i 



t ft ( 



o 



A 



o 



□, da 



1.0 1.5 2.0 

Log(3.6 jum Flux Density (/-tJy)) 



2.5 



FIG. 6. — The upper left quadrant of this plot indicates the infrared selection area for highly obscured, z > 2 type 2 AGN, defined by Marti'nez-Sansigre 
et al. (2005). We plot sources with > 5cr detections at 3.6 fim and 24 fim. Black circles surround all spectroscopically confirmed z > 2 sources. This selection 
quadrant successfully identifies all three SEXSI-confirmed z > 2 narrow-lined quasars. The ten spectroscopically unidentified SEXSI sources (no ID) in the 
selection quadrant are faint in the optical (typically R > 24) and have X-ray properties that are consistent with the notion that they are indeed high-redshift, 
obscured quasars. 



